Superfluid and Fermi liquid phases of Bose-Fermi mixtures in optical lattices 
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We describe interacting mixtures of ultracold bosonic and fermionic atoms in liarmonically con- 
fined optical lattices. For a suitable choice of parameters we study the emergence of superfluid 
and Fermi liquid (non-insulating) regions out of Bose-Mott and Fermi-band insulators, due to finite 
Boson and Fermion hopping. We obtain the shell structure for the system and show that angular 
momentum can be transferred to the non-insulating regions from Laguerre-Gaussian beams, which 
combined with Bragg spectroscopy can reveal all superfluid and Fermi liquid shells. 
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Fermi and Bose degenerate quantum gases and liquids 
are amazing systems, which have revealed individually 
several macroscopic quantum phenomena. For instance, 
superfluidity is known to exist in neutral liquids such as 
''He (boson) and in ''He (fermion), as well as in a vari- 
ety of electronic materials studied in standard condensed 
matter physics. The role of quantum statistics and inter- 
actions is of fundamental importance to understand the 
phases emerging from purely bosonic or purely fermionic 
systems, and a substantial amount of understanding of 
these individual Bose or Fermi systems can be found 
in the atomic and condensed matter physics literature. 
However, new frontiers can be explored when mixtures 
of bosons and fermions are produced in harmonic traps 
or optical lattices. An important example of the rich- 
ness of quantum degenerate Bose-Fermi mixtures was re- 
vealed in standard condensed matter systems, where for 
fixed ^He density and increasing amounts of ^He, the 
critical temperature for superfluidity is reduced, and be- 
low a tricritical point phase separation appears [H. In 
this Bose-Fermi mixture of standard condensed matter 
physics, essentially the only control parameter is the ra- 
tio between the densities of ''He and "^He. 

In atomic physics, a spectacular degree of control has 
been achieved in Bose-Fermi mixtures, where not only 
the ratio between densities of bosons and fermions can be 
adjusted, but also the interactions between fermions and 
bosons can be controlled through the use of Fcshbach res- 
onances , as demonstrated in mixtures of harmonically 
trapped Bose-Fermi polarized Fermion mixtures of '"^K 
and *^Rb. Furthermore, these same atoms have been suc- 
cesfully loaded into optical lattices Q and have produced 
a system that has no counterpart in standard condensed 
matter systems. By controlling the depths of optical lat- 
tices we can change not only the interactions between 
boson and fermions, but also their hopping from site to 
site, thus allowing the exploration of a very rich phase 
space, where supersolid and phase separated states have 
been suggested The list of Bose-Fermi mixtures in 
atomic physics is growing, and include systems where 
the masses are close like ^Li and ^Li, ^^K and ^°K, or 



^^^Yb and -'^''Yb; or systems where the masses are quite 
different like ^Li and ^^K, "^U and ^Sr, ^Li and ^SNa, 
^Li and s'^Rb, and ^SNa, or and ^'^Rb. This sug- 
gests a wide possibility of regimes that can be reached by 
tuning interactions, density and geometry, which is not 
is possible in ordinary condensed matter physics. 

A few studies of quantum phases of Bose-Fermi mix- 
tures have focused on homogeneous three dimensional 
systems with 0, 0] , and without optical lattices Q . How- 
ever, only one effort focused on harmonically confined 
optical lattices Q. Most of the descriptions of Bose- 
Fermi mixtures in optical lattices have relied on numer- 
ical methods using either Gutzwiller projection @| or 
quantum Monte Carlo 0] techniques. In this paper, we 
present a fully analytical theory of boson and spin po- 
larized fermion mixtures in harmonically confined opti- 
cal lattices by using degenerate perturbation theory for 
finite hopping in conjunction with the local density ap- 
proximation. This work provides insight into the phase 
diagram of Bose-Fermi mixtures, and into the detection 
of superfluid and Fermi liquid shells at low temperatures. 

This paper we analyses in detail the regime where the 
hopping parameters of bosons and fermions are compa- 
rable and the repulsion between bosons and fermions 
is a substantial fraction of the boson-boson repulsion. 
In this case, the system presents regions of (I) coexist- 
ing Bose-Mott and Fermi-band insulator, (II) coexisting 
Bose-Mott insulator and Fermi liquid, (HI) Bose-Mott 
insulator, and (IV) Bose superfluid, as shown in Fig. 1. 
We compute analytically the boundaries between various 
phases, and obtain the spatially dependent boson and 
fermions filling fractions in each region. Although one 
can envisage other situations where, for example, one 
can have coexistence of superfluid and Fermi band insu- 
lator, we confine our discussion to the situation above 
for the sake of simplicity. Finally, we propose a detection 
method of the shell strucuture of Bose-Fermi mixtures 
by using Lag uerre-Gaussian beams 0] and Bragg spec- 
troscopy [10[, where angular momentum is transferred 
only to regions with extended states such as the super- 
fluid and Fermi-liquid shells. 



2 



To describe Bose-Fermi mixtures in harmonically con- 
fined square (2D) or cubic (3D) optical lattices we start 
with the Hamiltonian 

H = /-Cb + A> - ^^s(r)7iB(r) - ^^i7(r)ni7(r) 



U 



BB 



^??.B(r) [ns(r) - 1] + C/BF^n-B(r)"-i^(r), 



where Kb = -ts I](r,r') and Kp = 

^^J" E(r r') /r/r' ^re the kinetic energies of boson 
and fermions with nearest-neighbor hoppings Ib and tp, 
and &J and fl are the bosonic and fcrmionic creation 
operators at site r. Here, the lattice sites for bosons and 
fermions are assumed to be the same, but the hopping 
parameters can be different. The number operators are 
nB(r) = h\hr and nF(r) = /r/r, and the corresponding 
local chemical potentials are ij,f{'^) = fJ-F — ^F(r) and 
A^s(r) = ^J'B — ^^B(r), where Vpir) = flpir/a)^ /2 and 
Vb{j^) = ^b{t / aY /2 are the harmonically confining 
potentials and and are the chemical potentials 
for fermions and bosons. The origin of the lattice 
with spacing a is chosen to be at the minimum of the 
harmonically confining potential. The terms containing 
Ubb {Ubf) represent the boson-boson (boson-fcrmion) 
interaction. 

When tB = tp = 0, the Hamiltonian H = Hq is 
a sum of single-site contributions, and the eigenstates 
are tensor products of number states with state vec- 
tors IV^) = \nBfi,nB,i,- ■ ■)\nF,o,nF^i,- ■ ■), with hb^t = 
0,1,2,... and np,. = 0,1 representing the occupation 
number of bosons and fermions at site r, respectively. 
At site r the local energy is £'„^^„^(r) ~ UBBnB{nB — 
1) /2 + UBFnBnp — l-i'B{i^)nB — fJ'F{i^)np. For the ground 
state wavefunction the number of bosons at site r is deter- 
mined by max{0, [{Ubb + f-iB{r))/2UBB\) if £'nB,o(r) < 
-E'„s,i(r) and max{0, [{Ubb + fJ'B{r)-UBF)/2UBB\) oth- 
erwise. Similarly, the number of fermions per site is zero 
if Engfi{r) < EnB,i{r) and one otherwise. The sym- 
bol [-J is the floor function. In the ground state solu- 
tion shells (nB^np) with hb bosons and np fermions 
are formed by those lattice sites r for which the lo- 
cal energy is the same. For our harmonic traps these 
shells are nearly spherically symmetric. The bound- 
ary between shells (uB^np) and {nB + l,np) is de- 
termined by E„^^„p{r) = ij;„^+i^„^ (r), leading to the 
radius R.B,nB.nF = ay/fln^^np/^B, where ft„s.np = 
2{fiB — ubUbb — upUbf)- Similarly, the bound- 
ary between shells with occupation numbers (nB,0) 
and (nB,l) is determined by equating the local ener- 
gies £'„^_„j^(r) and £'„g_„^_|_i(r), leading to the radius 
Rp.uB = a\/2(/ZF — nBUBF) I^F- We consider the num- 
ber of particles to be sufficiently large such that the radii 
of the boundaries are much larger than a. 

Next, we begin our discussion of finite hoppings by 
taking first Ib ^ 0, with tp — 0. The Bose su- 



perfluid region emerges due to kinetic fiuctuations at 
the boundaries between the {nB,np) and (n^ -I- l,np) 
shells. At this boundary the local energy i?„^+i „y(r) 
is degenerate with EnB,nF{i^)- To describe the emer- 
gence of superfluid regions, we introduce the order pa- 
rameter for superfluidity tpB.j via the transformation 

blbj — *■ i>*B^ihj + bji^Bj — tpB,i'^B.], and then for ana- 
lytical convenience make the continuum approximation 
ipir + a) = ip{r) + a.j9,V^(r) + {l/2)a^a.ididjil}{r). 

In the limit of Ubb ^ tB, we can restrict our 
Hilbert space to the number basis states \nB,np) and 
\nB + 1, np), as any contribution from other basis states 
to the local energy is of order t^/UBB- The hopping term 
tB affects the energies Ens+i,np{r:) and £;„g^„^(r) by 
removing their degeneracy, thus creating finite- width su- 
perfluid regions between shells (ns-l-l, np) and {nB,np). 
The effective local Hamiltonian then becomes, 



h: 



cff 



^n^ ,np(r) + A(r) -Vns + lA(r) 
-VnB + lA*(r) £;„^+i,„^(r)+A(r) 



(1) 

where A(r) = i(A(r)V'*(r) + cc) and A(r) = tB{zi'{r) + 
a^V^'0(r)). Here, z is the coordination number which 
depends on the lattice dimension d. 

The eigenvalues of Eq. ([1]) are given by, 

E±{v) = E,{v) ± ^[Ea{v)f + {nB + l)\m\\ 

where i;,(r) = [i;„^+i,„^) (r) + £;„^,„^ (r)] /2 A(r) 
is proportional to the sum of the diagonal terms, and 

Ed{v) = [EnB + l^npi-r)- EnB,nF{-r)\/2 = {nBUBB + 

npUBF~ M-b(i")/2 is proportional to their difference. No- 
tice that E-(y) is the lowest local energy leading to the 
total ground state energy E = -p- / (iri?_(r). 

The order parameter equation (OPE) is determined by 
minimization of E with respect to '0*(r) leading to 



A(r) 



(nB + l)tB(^ + a'V')A(r) 
2y'|i?4r)| V (ns + 1) |A(r)r 



0. 



(2) 



Notice that the OPE is not of the Gross-Pitaeviskii (GP) 
type, since the superfluid regions emerge from local fluc- 
tuations between neighboring Mott shells. Ignoring the 
spatial derivatives of ij) in Eq. [5] leads to the spatially 
dependent order parameter 



2 nB + l [nBUBB + npUBF - HB{'r)Y 



izHj^inB + 1) 



(3) 



Since \ip{r)\ > 0, hence \nBUBB + npUBF — ^J■B,r\ < 
(riB + l)ztB; and the inner i?„g „p^_ and outer i?„g „j^ .|_ 
radii for the superfluid shell between the {riB,np) and 
(ns -I- l,np) Mott regions are obtained by setting 
\'4'(y)\ =0 leading to 



R 



'i± 



2ztB{nB + 1) 



d2 

B,nB,nF 
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This relation shows expUcitly that ts sphts the spatial 
degeneracy of the [nB,np) and [ub + Ij't.f) insulating 
shells at r Rc,nB,nF or = nsUBB + upUbf 

by introducing a superfluid region of width Ai?„g „p = 
RnB,nF,+ — RnB,np,-- (See Fig. [T] for characteristic 
widths) . 

In addition, the local bosonic filling fraction 

dE^{v) ^1 ubUbb + tifUbf - 

2' 



ns(r) = - 



dfiB ' 2 2ztB{nB + 1) 

in the same region interpolates between hb + 1 for 
r ^ RnB,np- and n for r > RnBMF,+ , while the 
chemical potential hb is fixed by the total number of 
bosons Nb = J drn(r). The local bosonic compressibil- 
ity KB{r) = dnB{r)/dnB = ^/2ztB{nB + 1) of the super- 
fluid shells is non-zero, in contrast to the incompressible 
(kb = 0) {nB,np) and [ub + 1,?^_f) insulating shells for 
r < RnB.riF- and r > Rns,np,+ , respectively. 
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FIG. 1; (color online) (a) Shell structure of Bose-Fermi mix- 
tures in harmonically confined optical lattices showing a coex- 
isting Bose-Mott and Fermi-band insulator region at the cen- 
ter (n_B = IjfiF = 1), a coexisting Bose-Mott insulator and 
Fermi liquid region (in blue), a Bose-Mott insulator region 
(ns = IjUf = 0), and Bose-Superfluid region at the edge 
(in red), (b) filling factors for fermions shown as the solid 
(dark blue) curve and for bosons shown at the dashed (red) 
curve. The solid parabolic curve (light blue) shows the or- 
der parameter in the superfiuid region. The parameters used 
are tp = ts = 0.0325Ubb, Ubf = Q.IUbb, Ms = 0.8l7ss, 
/If = QAUbb and Q,f = Os = 8 x 1Q~^Ubb, which are rep- 
resentative of Bose-Fermi mixtures with nearly the same mass 
such as "^Li and ^Li, =^^K and "'"K, or ^^^Yb and ^'^^Yb. For 
the parameters chosen, the widths of the superfiuid and FL 
shells are several times larger than the lattice spacing a. 

Now, we consider finite tp- In order to have a tractable 
theory we assume that the shell boundaries of the bosons 
and fermions are well separated. This allows us to inves- 
tigate the Fermi liquid near the shell boundary of the 
fermions in the presence of a Bose-Mott insulator with 
Ub bosons per site. Furthermore, if we assume that the 
local density of the Fermi gas is smoothly varying then 
the local number of fermions is 



np{Y) 



(4) 



where f[x] is the Fermi function at temperature T, and 
2?(e) = ^jj. S{e — Ek) is the density of fermion states with 



energy dispersion ek = —'2tpJ2i cos{kia). The effective 
chemical potential /i|?(r) = ^J,p{r) — ubUbf accounts 
for the effect of the bosons. The band minimum and 



maximum of ek are Cn 



-2dtp and 



2dtf 



respectively. Thus, the Fermi liquid region is limited 
by the boundaries emin < A*-F(r) 

Rf,± = RF,n 



< Cmax, leading to 
^1 ± 2ztpa? /Q^pR^p for the inner 
and outer radius of the FL shell. The width of the 

FL region is ARp ^B — Rf,+ — Rf,-- (See Fig.[l]for char- 
acteristic widths) . The isothermal compressibility of the 
FL region is hf{y) = dnp{r)/diJ.F, which leads at zero 
temperature to ^^(r) = outside the FL shell, indicating 
the presence of insulating regions and ^^(r) = ©[//^(r)] 
inside the FL shell, indicating the presence of conducting 
regions. The superfluid and FL shells for finite is and 
tp, and their density profiles are shown in Fig. [T]for the 
two-dimensional case. 

Next, we propose an experiment to detect superfluid 
and Fermi liquid shells in Bose-Fermi mixtures using a 
combination of Gaussian and Laguerre-Gaussian beams 
followed by Bragg spectroscopy. To illustrate the idea, we 
discuss the simpler case of a nearly two-dimensional con- 
flguration, where the harmonic trap is very tight along 
the z-direction, loose along the x- and y- directions. Upon 
application of Gaussian and Laguerre-Gaussian beams 
along the z-direction, only angular momentum is trans- 
ferred to the atoms in the conducting phases (superfluid 
or Fermi liquid) , imposing a rotating current with a well 
defined velocity profile, while the insulating regions do 
not absorb angular momentum due to their large gap in 
the excitation spectrum. 




FIG. 2: (Color online) Schematic plot for the detection of 
outer (red) superfiuid and inner (blue) Fermi liquid shells us- 
ing Bragg spectroscopy. The angles 9i and 62 indicate the 
locations of strongest momentum transfer from the Bragg 
beams (large green arrows) to the rotating superfluid and 
Fermi liquid shells of radii Ri and i?2. The gray arrows indi- 
cate the sense of rotation of the conducting shells. 

To probe the rotating superfluid and Fermi liquid 
phases we propose the use of two counter-propagating 
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Bragg beams applied in the xy plane along the x direc- 
tion, as indicated in Fig. [2l 

The Bragg beams transfer a net linear momentum 
h{k + fc')x to the atoms of mass m which satisfy the 
energy conservation condition 

h(Lj - u}') =ef - ei- v^hlk + k') + M^LtAi_^ (5) 

2m 

where Vx is the component of the velocity v(r) — p(r)/m 
along the x direction and and Cf are the experimen- 
tally accessible energies of the initial and final internal 
states atom. For an atom carrying one unit of angu- 
lar momentum, the velocity is = hO/mr. Therefore, 
within a conducting shell with radius r ^ R atoms get 
a linear momentum kick of H{k + k')'k when the velocity 
Vx = hsinO/mR satisfies the condition given in Eq. [5l 
This leads to two Bragg angles 9 = — sm~^ {mRvx / h) , 
and TT — 6 for each conducting shell. As can be seen in 
Fig. [21 the Bragg angles are 9i and tt — for the outer 
superfluid shell labelled by Ri , and arc 62 and tt — 62 for 
the fermi liquid shell labelled by i?2- Once these atoms 
are kicked out of the conducting shells, they form two 
small expanding clouds, which can be detected by direct 
absorption imaging. 

Next, we discuss the time scales over which the ro- 
tation in the conducting regions persist and can be de- 
tected experimentally. In the case of the superfluid region 
we use the Landau criterion to show that the velocity 
imposed to the superfluid through the angular momen- 
tum transfer is much smaller than the local sound ve- 
locity c(r) — y/ps{r)/K, where ps(r) = 2t_Ba^|V'(r)P is 
the local superfluid density, and k is the compressibility. 
Thus, c(r) ~ 2y/ {ub + l)zta\ip{r)\ vanishes at the insu- 
lator boundaries where |'0(r)| = 0, and only close to the 
edge of the superfluid regions the local rotational speed 
v{r) = h/mr exceeds c(r), which means that essentially 
all the superfluid region can be detected and the angular 
momentum transferred docs not decay over time scales 
of at least seconds, limited by the lifetime of the trapped 
system. 

In the case of the Fermi liquid region, the time scale 
over which the flow of the fermions persist in presence 
of the Bosc-Mott insulator background can be calculated 
from the imaginary part of the fermionic self-energy 

^F{k) = UlpT^ J2 GB{qi)GB{q2)GF{k + qi - 92), 

91,92 

where k = (k, zw) and qi = {cii,ivi), with uj (vi) are 
fermionic (bosonic) Matsubara frequencies and T is tem- 
perature. The bare inverse bosonic propagator in the 
Bose-Mott phase is 



1 



ns + 1 



Ub 



ihuj - 



where Ei{t) = [ub 
2)Ubb - MB(r) and Eq 



£;i(r) ihuj~E2{r) 

1)C^BS - £'2(r) = {ub - 

cos(A:^a). The bare 



inverse fermionic propagator in the Fermi liquid phase 
is Gp^{k, r) = ihuj — epi^, r), where eF(k, r) = eF(k) — 
fip{r). For Ub = 1 and T ~ 0, the imaginary part of the 
fermionic self-energy is 



ImI]F(fc,r) ^ -ttUIp [F{huj) + F{~huj)] , 



(6) 



where F{nuj) = Q {nw)e {- tiuj + ^B{r))V{huj + urir)), Q 
is the Heaviside step function, and I?(e) is the Fermion 
density of states. For a two-dimensional Fermi liquid 
shell, there is a Van Hove singularity in I?(e) at half flU- 
ing. The expression in Eq. ^ is independent of mo- 
mentum, since the dominant excitations in the Bose- 
Mott region are number-conserving and low-momentum 
particle- hole excitations, but strongly dependent on posi- 
tion through /iB(r) and /i_F(r), leading to a characteristic 
decay time ~ — /i/lml](/s, r). For the parameters used 
in Fig. [2 (with Usp/h = 1 kHz) the time scale for the 
persistence of the flow near the edges (away from the Van 
Hove singularity) is ~ 13 ms. However, near the center 
of the Fermi liquid region (close to the Van Hove singu- 
larity) Tr is extremely short, indicating that it is much 
easier to detect fermions at the edge than at the center 
of Fermi liquid shells. 

We have discussed the phase diagram of Bosc-Fcrmi 
mixtures in harmonically confined optical lattices in 
the regime where the hopping parameters of bosons 
and fermions are comparable and the repulsion between 
bosons and fermions is a substantial fraction of the 
boson-boson repulsion. We showed that the system ex- 
hibits regions of (I) coexisting Bose-Mott and Fermi- 
band insulator, (II) coexisting Bose-Mott insulator and 
Fermi liquid, (HI) Bose-Mott insulator, and (IV) Bosc- 
superfluid. We have calculated analytically the bound- 
aries between these phases and obtained the spatially de- 
pendent filling fraction for each region. Finally, we pro- 
posed a detection method of the superfluid and Fermi 
liquid shells of Bose-Fcrmi mixtures by using Gaussian 
and Lagucrre-Gaussian beams followed by Bragg spec- 
troscopy. 
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